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Motivated by recent experiments on quasi-lD vanadium oxides, we study quantum phase transi- 
tions in a one-dimensional spin-orbital model describing a Haldane chain and a classical Ising chain 
locally coupled by the relativistic spin-orbit interaction. By employing a field-theoretical approach, 
we analyze the topology of the ground-state phase diagram and identify the nature of the phase 
transitions. In the strong coupling limit, a long-range Neel order of entangled spin and orbital 
angular momentum appears in the ground state. We find that, depending on the relative scales 
of the spin and orbital gaps, the linear chain follows two distinct routes to reach the Neel state. 
First, when the orbital exchange is the dominating energy scale, a two-stage ordering takes place 
in which the magnetic transition is followed by melting of the orbital Ising order; both transitions 
belong to the two-dimensional Ising universality class. In the opposite limit, the low-energy orbital 
modes undergo a continuous reordering transition which represents a line of Gaussian critical points. 
On this line the orbital degrees of freedom form a Tomonaga-Luttinger liquid. We argue that the 
emergence of the Gaussian criticality results from merging of the two Ising transitions in the strong 
hybridization region where the characteristic spin and orbital energy scales become comparable. 
Finally, we show that, due to the spin-orbit coupling, an external magnetic field acting on the spins 
can induce an orbital Ising transition. 



I. INTRODUCTION 

Over the past decades, one-dimensional spin-orbital 
models have been a subject of intensive theoretical stud- 
ies. The interest is to a large extent motivated by exper- 
imental discovery of unusual magnetic properties in vari- 
ous quasi-one-dimensional Mott insulators The inter- 
dependence of spin and orbital degrees of freedom is usu- 
ally described by the so-called Kugel-Khomskii Hamil- 
tonian in which the effective spin exchange constant de- 
pends on the orbital configuration and vice versa^i^ An- 
other mechanism of coupling spin and orbital degrees of 
freedom is the on-site relativistic spin-orbit (SO) interac- 
tion A L-S, where L is the orbital angular momentum and 
A is the coupling constant. In compounds with quenched 
orbital degrees of freedom, the presence of the SO term 
usually leads to the single-ion spin anisotropy DS^ where 
D ~ A 2 /A and A denotes the energy scale of the crystal 
field which lifts the degenerate orbital states. 

For systems with residual orbital degeneracy, on the 
other hand, the effect of the SO term is much less 
explored compared with the Kugcl-Khomskii-type cou- 
pling. Due to the directional dependence of the orbital 
wave functions, the SU(2) symmetry of the Heisenberg 
spin exchange is expected to be broken in the presence of 
the SO interaction. The resultant spin anisotropy is likely 
to induce a long-range magnetic order in the spin sector. 
A more intriguing question is what happens to the orbital 
sector. To answer this question, one needs to consider 
the details of the interplay between the orbital exchange 
and the SO coupling. Here we consider the simplest 



case of a two-fold orbital degeneracy per site. Specifi- 
cally, the two degenerate states could be the d yz and d zx 
orbitals in a tetragonal crystal field observed in several 
transition-metal compounds. We introduce pseudospin- 
1/2 operators r a (a = x,y,z) to describe the doublet 
orbital degrees of freedom assuming that t z = ±1 corre- 
spond to the states \yz) and i\zx), respectively. Alterna- 
tively, one can also realize the double orbital degeneracy 
in the Mott-insulating phase of a ID fermionic optical 
lattice where the eigenvectors of r z refers to p x and p y 
orbitals in an anisotropic potential££ Restricted to this 
doublet space, the orbital angular momentum operator 
L = (0, 0, t x ). This can be easily seen by noting that the 
eigenstates of t x carry an angular momentum (L z ) = ±1. 

The exchange interaction between localized orbital de- 
grees of freedom is characterized by its highly direc- 
tional dependence: the interaction energy only depends 
on whether the relevant orbital is occupied for bonds of 
a given orientation. This is particularly true for inter- 
actions dominated by direct exchange mechanism. De- 
noting the relevant orbital projectors on a given bond as 
P = (l+r^/2, where is an appropriate pseudospin- 
1/2 operator (r' 9 being a Pauli matrix), the orbital inter- 



action is thus described by an Ising- type term t^Tj. 
well studied orbital compass model and Kitaev model 
both belong to this category^iS The quantum nature of 
these models comes from the fact that different operators 
t' 3 , which do not commute with each other, are used for 
bonds of different types. To avoid unnecessary compli- 
cations coming from the details of orbital interactions, 
we assume that there is only one type of bond in our ID 
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system and the orbital interaction is thus governed by a 
classical Ising Hamiltonian. 

We incorporate these features into the following toy 
model of spin-orbital chain ( J s , J T > 0): 

H = H s + H t + H St (1) 

= Js 2J S " ' S »+l + ^ T n T n+l + ^ T n^n- 

n n n 

Motivated by the recent experimental characterizations 
of quasi- ID vanadium oxides,— ~— here we focus on the 
case of quantum spin with length S = 1. The above 
model thus describes a Haldane chain locally coupled to a 
classical Ising chain by the SO interaction Hs T - The role 
of the A-term is two-fold: firstly it introduces anisotropy 
to the spin-1 subsystem, and secondly it endows quantum 
dynamics to the otherwise classical Ising chain. 

Before turning to a detailed study of the phase dia- 
gram of model (p}, we first discuss its connections to 
real compounds. As mentioned above, the interest in 
the toy model is partly motivated by the recent experi- 
mental progress on vanadium oxides which include spinel 
ZnV 2 0^i-il and quasi-lD CaV 2 4 3^ In both types 
of vanadates, the two d electrons of V 3+ ions have a 
spin S = 1 in accordance with Hund's rule. In the low- 
temperature phase of both vanadates, the vanadium site 
embedded in a flattened VOg octahedron has a tetragonal 
symmetry. This tetragonal crystal field splits the degen- 
erate t2 g triplet into a singlet and a doublet. As one of 
the two d electrons occupies the lower-energy d xy state, 
a double orbital degeneracy arises as the second electron 
could occupy either d zx or d yz orbitals. The fact that 
the d xy orbital is occupied everywhere also contributes to 
the formation of weakly coupled quasi-lD spin-1 chains in 
these compounds^ On the other hand, the details of the 
orbital exchange depends on the geometry of the lattice 
and in the case of vanadium spinel the orbital interaction 
is of three-dimensional nature. The Ising orbital Hamil- 
tonian in Eq. (fTJ) thus should be regarded as an effective 
interaction in the mean-field sense. Nonetheless, the toy 
model provides a first step towards understanding the es- 
sential physics introduced by the SO coupling. Moreover, 
many conclusions of this paper can be applied to the case 
of quasi-lD compound CaV204 where the vanadium ions 
form a zigzag chain. 

It is instructive to first establish regions of stable mas- 
sive phases. In the decoupling limit, A — > 0, our model 
describes two gapped systems: a quantum spin-1 Heiscn- 
berg chain and a classical orbital Ising chain. The ground 
state of the spin sector is a disordered quantum spin 
liquid with a finite spectral gap2i A5, whereas the or- 
bital ground state is characterized by a classical Neel 
order along the chain: (r*) = (— l) n rj z . Quantum ef- 
fects in the orbital sector induced by the SO coupling 
play a minor role. Obviously, just because of being 
gapped, both the spin-liquid phase and the orbital or- 
dered state are stable as long as A remains small. Con- 
sider now the opposite limit, A 3> Js, J T - In the ze- 
roth order approximation, the model is dominated by the 



single-ion term Hs T whose doubly degenerate eigenstates 
|±) = \S Z = ±1) <£> \t x = ±1) represent locally entangled 
spin and orbital degrees of freedom. Switching on small 
Js and J T leads to a staggered ordering of the |+) and 
|— ) states along the chain. Physically, the large-A ground 
state can be viewed as a simultaneous Neel ordering of 
spin and orbital angular momentum characterized by or- 
der parameters ( and rj x such that (S%) = (— 1)"C an d 
(LfJ = (t x ) = ( — \) n rj x . The Ising order parameter rf 
vanishes identically in this phase. 

These observations naturally lead to the following 
questions. How is the magnetically ordered Neel state 
at large A connected to the disordered Haldane phase 
as A — > ? What is the scenario for the orbital re- 
orientation transition rf if , which is of essentially 
quantum nature ? In this paper we employ the field- 
theoretical approach to address these questions. We first 
note that the one-dimensional model (|TJ) is not exactly 
integrable. As a consequence, the regime of strong hy- 
bridization of the spin and orbital excitations, which is 
the case when J T , Js and A are all of the same order, 
stays beyond the reach of approximate analytical meth- 
ods. We thus will be mainly dealing with limiting cases 
J T ^> J s and J T <^ Js, in which one can integrate out 
the "fast" variables to obtain an effective action for the 
"slow" modes. Following this approach, we establish the 
topology and main features of the ground-state phase 
diagram in the accessible parts of the parameter space 
of the model. We were able to unambiguously identify 
the universality classes of quantum criticalitics separat- 
ing different massive phases. Using plausible arguments 
we comment on some features of the model in the regime 
of strong spin-orbital hybridization. 

We demonstrate that the aforementioned reorientation 
transition rf —¥ rf can be realized in one of two possi- 
ble ways. In the limit of large J T , we find a sequence 
of two quantum Ising transitions and an intermediate 
massive phase, sandwiched between these critical lines, 
in which both rf and rf are nonzero. This is consis- 
tent with the recent findings^ based on DMRG calcu- 
lations and some analytical estimations. In the oppo- 
site limit, when the Haldane gap Ag is the largest en- 
ergy scale, integrating out the spin excitations yields an 
effective lowest-energy action for the orbital degrees of 
freedom, which shows that the rf — > rf crossover takes 
place as a single Gaussian quantum criticality. At this 
critical point, the orbital degrees of freedom display an 
extremely quantum behaviour: they arc gapless and form 
a Tomonaga-Luttinger liquid. This is the main result of 
this paper. We bring about arguments suggesting that 
the emergence of the Gaussian critical line is the result 
of merging of the two Ising criticalities in the region of 
strong spin-orbital hybridization. 

Any field-theoretical treatment of the model (JTJ) must 
be based on a properly chosen contiuum description of 
the spin-1 antiferromagnetic Hciscnbcrg chain. Its prop- 
erties have been thoroughly studied, both analytically 
and numerically (see for a recent review Ref. 1231 ). In 
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what follows, the spin sector of the model (QJ will be 
treated within the 0(3)-symmetric Majorana field the- 
ory, proposed by Tsvelik^ 
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Here £,%l(x) is a degenerate triplet of real (Majorana) 
Fermi fields with a mass m, the indices R and L label 
the chirality of the particles, and 
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is a weak four-fcrmion interaction which can be treated 
perturbatively. The continuum theory ([2]) adequately de- 
scribes the low-energy properties of the generalized spin-1 
bilinear-biquadratic chain 



Hs^Hs 



'n+l 



(3) 



in the vicinity of the critical point j3 = 1.— This quantum 
criticality belongs to the universality class of the SU(2)2 
Wess-Zumino-Novikov-Witten (WZNW) model with cen- 
tral charge c = 3/2. 

At small deviations from criticality the Majorana mass 
m ~ Js\P ~ 1| determines the magnitude of the triplet 
gap, As = \m\ <C Js- The theory of a massive 
triplet of Majorana fermions is equivalent to a system 
of three degenerate noncritical 2D Ising models, with 
m ~ (T — T c )/T c . This is one of the most appealing fea- 
tures of the theory because the most strongly fluctuating 
physical fields of the 5 = 1 chain, namely the staggered 
magnetization and dimerization operators, have a simple 
local representation in terms of the Ising order and dis- 
order parameters i 24 i 26 ' 27 It is this fact that greatly sim- 
plifies the analysis of the spin-orbital model ([T]). While 
the correspondence between the models ([2]) and (|3|) is 
well justified at \/3 — 1| <C 1, it is believed that the Ma- 
jorana model © captures generic properties of the Hal- 
dane spin- liquid phase of the spin-1 chain, even though at 
large deviations from criticality — 1| ~ 1, A5 ~ J s ) 
all parameters of the model should be treated as phe- 
nomcnological. 

The remainder of the paper is organized as follows. We 
start our discussion with Sec. [TT| which contains a brief 
summary of known facts about the Majorana modcl 2 ^ 
that will be used in the rest of the paper. In Sec. IIIII 
we consider the limit J T /As S> 1 and by integrating 
out the 'fast' orbital modes, show that on increasing the 
SO coupling A the system undergoes a sequence of two 
consecutive quantum Ising transitions in the spin and or- 
bital sectors, respectively. In section IIVI we analyze the 
opposite limiting case, J T /As <C 1, and, by integrating 
over the 'fast' spin modes, show that there exists a sin- 
gle Gaussian transition in the orbital sector accompanied 
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FIG. 1: Schematic phase diagram of the model on the (xs, 
x T )-plane, where xs = As/A and x T — J T /A. 



by a Neel ordering of the spins. We then conjecture on 
the topology of the ground-state phase diagram of the 
model. In Sec. [V] we show that spin-orbital hybridiza- 
tion effects near the orbital Gaussian transition lead to 
the appearance of a non-zero spectral weight well below 
the Haldane gap which can be detected by inelastic neu- 
tron scattering experiments and NMR measurements. In 
Sec. I VII wc comment on the role of an external magnetic 
field. We show that, through the SO interaction, a suf- 
ficiently strong magnetic field affects the orbital degrees 
of freedom and can lead to a quantum Ising transition 
in the orbital sector. Sec. IVIII contains a summary of 
the obtained results and conclusions. The paper has two 
appendices containing certain technical details. 



II. SOME FACTS ABOUT MAJORANA 
THEORY OF SPIN-1 CHAIN 

In this Section, we provide some details about the 
0(3)-symmetric Majorana field theory^ Eq. @, which 
represents the continuum limit of the biquadratic spin-1 
model © at \j3 - 1| < 1. 

In the continuum description, the local spin density of 
the spin model §3§ has contributions from the low-energy 
modes centered in momentum space at q = and q = it: 



S(x) = I R (x) + I L (x) + (-iy/ a °N(x) 



(4) 



The smooth part of the local magnetization, I = Ir + 
II, is a sum of the level-2 chiral vector currents. The 
SU(2)2 Kac-Moody algebra of these currents is faithfully 
reproduced in terms of a triplet of massless Majorana 
fields^ I = (£\£ 2 ,£ 3 ): 



(5) 
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This fact is not surprising because, as already men- 
tioned, the central charge of the SU(2) 2 WZNW theory 
is c = 3/2, whereas that of the theory of a massless Ma- 
jorana fermion (equivalently, critical 2D Ising model) is 
c = 1/2. At small deviations from criticality (|/3 — 1| -C 1) 
the fermions acquire a mass. Strongly fluctuating fields of 
the spin-1 chain, the staggered magnetization N(x) and 
dimerization operator e(x) = (— 1)™S„ • S n +i, arc nonlo- 
cal in terms of the Majorana fields but admit a simple 
representation in terms of the order, a, and disorder, //, 
operators of the related noncritical Ising models: 



N ~ (1/a) (<7i/i 2 M3i ^102^3, M1M203), 
e ~ (l/a)cri (T2(T3, 



(0) 



where a ~ a is a short-distance cutoff of the contin- 
uum theory These expressions together with their duals 
(i.e. their counterparts obtained by the duality trans- 
formation in all Ising copies, a a O fXa) determine the 
vector and scalar parts of the WZNW 2x2 matrix field 
g which is a primary scalar field with sca ling dimension 
3/8. It has been demonstrated in Ref. [28| that using 
the representation ^ and the short-distance operator 
product expansions for the Ising fields, one correctly re- 
produces all fusion rules of the SU(2)2 WZNW model. 
An equivalent way to make sure that this is indeed the 
case is to consider the four-Majorana representation of 
the weakly coupled spin- 1/2 Heisenberg ladde r 26 ' 27 and 
take the limit of a infinite singlet Majorana mass to map 
the low-energy sector of the model on the 0(3) theory 

In the spin-liquid phase of the spin chain (|3]), which 
is the case j3 < 1, the Majorana mass m is positive, 
implying that the degenerate triplet of 2D Ising models is 
in a disordered phase: (a a ) = 0, (fi a ) 7^ (a = 1,2,3). In 
particular, this implies that the 0(3) symmetry remains 
unbroken, (N) = 0, and the ground state of the system 
is not spontaneously dimerized, (e) =0. 

The representation ([6]) proves to be very useful for cal- 
culating the dynamical spin correlation functions because 
the asymptotics of the Ising correlators (ct(x, r)cr(0, 0)) 
and (fJ.(x, t)/i(0, 0)) are well known both at criticality and 
in a noncritical regime. In the disordered phase (m > 0), 
the leading asymptotics of the Ising correlators are: 



(MrMO)} 
((7(r)«7(0)) 



K£s) 1/4 

(a/£s) V4 y/&/r , 



0( e -2r/?s 
-r/ts 



(7) 



where £s = v/m is the correlation length, and r — 
\f x 2 + v 2 t 2 . (By duality, in the ordered phase (m < 0) 
the asymptotics of the correlators in ([7]) must be inter- 
changed.) Correspondingly, the dynamical correlation 
function 



<N(r)N(0)) 



-r/£s 



(8) 



Its Fourier transform at q ~ 7r and small uj describes a 
coherent excitation - a triplet magnon with the mass gap 



3to x((7,w) ~ t — -5 (ui — y/(q — n) 2 v 2 + m 2 ) . (9) 

Since the single-ion anisotropy i? an is = DJ2 n (^n) 2 
lowers the original 0(3) symmetry down to 0(2) x Z2, 
one expects 2 ^ that in the continuum theory it will induce 
anisotropy in the Majorana masses 

mi = ma / m 3 , 

as well as in the coupling constants parametrizing the 
four-fcrmion interaction: 



Hi 



a^b 



This can be checked by using the correspondence ((31) 
and short-distance operator product expansions (OPE) 
for the physical fields. There will also appear anisotropy 
in the velocities, V\ = v- 2 7^ U3, but we will systematically 
neglect this effect. Thus, we have i? an is = J dx "Hanis, 
with 



= Da dx [l 3 {x)I 3 {x + a) + N 3 (x)N 3 (x + a)] , (10) 



where a ~ a is a short-distance cutoff of the continuum 
theory. Using ^ and keeping only the Lorentz invariant 
terms (i.e. neglecting renormalization of the velocities) 
we can replace (I 3 ) 2 by 2/| I 3 . To treat the second term 
in the r.h.s. of (jlOp . we need OPEs for the products of 
Ising operators^ 



<?{z, z)a(w, w) 
= ±f- 



V2 



/i(z, z)fi(w, w) 



\ 1/4 

) [l — ir\z — w\s(w, «;)] , (11) 



-31(1^)^* -•*<».•>]■(«> 

Here e = i£,Ft£,L is the energy density (mass bilinear) of 
the Ising model, z = vt + ix and w = vt' + ix' are 
two-dimensional complex coordinates, z and w are their 
conjugates. From the above OPEs it follows that 



N 3 {x)N 3 {x + a) = ifr/a) + - 

where C ~ 1 is a nonuniversal constant. As a result, 



H anis = -i 5m « + \ E St* (&&) , (13) 



0=1,2,3 



where 



Smi = 5rri2 = —Sm^ = —(irC)D 



(14) 
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are corrections to the singlc-fcrmion masses, and Sgi 2 = 
(2 — ir 2 C)Da 7 Sgi 3 = Sg 23 = ir 2 CDa are cou- 
pling constants of the induced interaction between the 
fermions. Smallness of the Majorana masses (|m|a/v <C 
1) implies that the additional mass renormalizations 
caused by the interaction in (|13[) are relatively small, 
m(Da/v) ln(v/\m\a) <C D, so that the main effect of 
the single-ion anisotropy is the additive rcnormalization 
of the fcrmionic masses, m a = m + 5m a , with 8m a given 
by Eq.CU. 

The cases D > and D < correspond to an easy- 
plane and easy-axis anisotropy, respectively. The spin 
anisotropy (|18j) induced by the spin-orbit coupling is of 
the easy-axis type. At D < the singlet Majorana 
fermion, £ 3 , is the lightest, m 3 < mi = m 2 . Increas- 
ing anisotropy drives the system towards an Ising crit- 
icality at D = -D*, where m 3 = 0. At D < -D* 
the system occurs in a new phase where the Ising dou- 
blet remains disordered while the singlet Ising system 
becomes ordered. It then immediately follows from the 
representation © that the new phase is characterized by 
a Neel long-range order with (N 3 ) ^ 0. Transverse spin 
fluctuations, as well as fluctuations of dimerization, are 
incoherent in this phase. 



III. TWO ISING TRANSITIONS IN THE 
A s <Jr LIMIT 

Now we turn to our model |T]). Let us consider the case 
when, in the absence of spin-orbit coupling, the orbital 
gap is the largest: J T ^> J s . The orbital pseudospins then 
represent the 'fast' subsystem and can be integrated out. 
Assuming that A <C J T , we treat the spin-orbit coupling 
perturbatively. In this case, the zero order Hamiltonian 
Ho = Hs + H T describes decoupled spin and orbital sys- 
tems, while the spin-orbit interaction Hs T denotes per- 
turbation. Defining the interaction representation for all 
operators according to A(t) 



rH 



°Ae~ THo (here r de- 
notes imaginary time), the interaction term in the Eu- 
clidian action is given by 



(15) 



The first nonvanishing correction to the effective action 
in the spin sector is of the second order in A: 



AS S 



-A 2 



nm 

(16) 



Averaging in the right-hand side of (fT6|) goes over config- 
urations of the classical Ising chain H T . The correlation 
function (t^(ti)t^(t 2 )) t is calculated in Appendix |A"1 It 
is spatially ultralocal (because there are no propagating 
excitations in the classical Ising model) and rapidly de- 
caying at the characteristic time ~ 1/Jt, which is much 



shorter than the spin correlation time ~ 1/Ao: 

«(nX(r a )) T = L cxp (-4J T \n - r 2 |) . (17) 

Passing to new variables, r = (n +r 2 )/2 and p = T\ — t 2 , 
and integrating over p yields a correction to the effec- 
tive spin action which has the form of a single-ion spin 
anisotropy. Thus in the second order in A, the spin 
Hamiltonian acquires an additional term 

-^rX [slf - (18) 



The anisotropy splits the Majorana triplet into a doublet 
K 1 !? 2 ) an d singlet (£ 3 ), with masses 



ttCX 2 



ttCX 2 



mi = n%2 = m + , m 3 = m — — , (19) 

where C ~ 1 is a nonuniversal positive constant. The 
anisotropy is of the easy-axis type, so that the singlet 
mode has a smaller mass gap. 

As long as all the masses m a remain positive, the sys- 
tem maintains the properties of an anisotropic Haldane's 
spin-liquid. The dynamical spin susceptibilities calcu- 
lated at small L) and q ~ n (see Sec. [TT| . 



mi 



j-<5 ^cj - y 1 {q- tt) 2 v 2 + mf^j 



3m 



indicate the existence of the S z = ±1 and S z — optical 
magnons with mass gaps mi and r/13, respectively. In- 
creasing the spin-orbital coupling leads eventually to an 
Ising criticality at A = A c i = I^J ' J T m/irC, where m 3 = 0. 
At m 3 < the system occurs in a long-range ordered Neel 
phase with staggered magnetization (S^) = (— l) n C(<\); 
in which the ^-symmetry of model (|18j) is spontaneously 
broken. Using the Ising-model representation ([5]) of the 
staggered magnetization of the spin-1 chain, we find that 
at < A — A c i <C A c i the order parameter £(A) follows a 
power-law increase: 



C(A) 



A — A c i 
Ad 



1/8 



(21) 



The transverse spin fluctuations become incoherent in 
this phase. The situation here is entirely similar to that 
in the spontaneously dimerized massive phase of a two- 
chain spin- 1/2 ladde r 27 i 30 , where the dimerization kinks 
make spin fluctuations incoherent. In the present case, 
the spontaneously broken Z 2 symmetry of the Neel phase 
leads to the existence of pairs of massive topological kinks 
contributing to a broad continuum with a threshold at 
ui = mi + |r7^3 1 (the details of calculation can be found in 
Refill): 

Sm X xx (q,u) (22) 
1 6{uj 2 - (q - tt)V - (mi + |m 3 |) 2 ) 



^mi\m 3 \ ^Juj 2 - {q - rr) 2 v 2 - (m 1 + |m 3 |) 5 
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In the Neel phase, the orbital sector acquires quantum 
dynamics because antiferromagnetic ordering of the spins 
generates an effective transverse magnetic field which 
transforms the classical Ising model H T to a quantum 
Ising chain. At A > A c i the spin-orbit term takes the 
form 

H S r = -hJ2(-^ + H' ST , (23) 

n 

where h = \((\) and H' Sj = -A£ n (S*- ac- 
counts for fluctuations. Since both the orbital and spin 
sectors are gapped, the main effect of this term is a 
renormalization of the mass gaps and group velocities. 
The transverse field h gives rise to quantum fluctuations 
which decrease the classical value of rf and, at the same 
time, lead to a staggered ordering of the orbital pseu- 
dospins in the transverse direction. Since the orbital sec- 
tor has a finite susceptibility with respect to a transverse 
staggered field, in the right vicinity of the critical point 
if follows the same power-law increase as £ but with a 
smaller amplitude: 




This result is in a good agreement with previously 
obtained numerical results for order parameters (See 
Fig. 4(a) in Ref.Sl). 

Performing an inhomogeneous 7r-rotation of the pseu- 
dospins around the y-axis, t%' z — > (— l) n T*' z , t% — > rjjf, 
we find that at A > A c i the effective model in the orbital 
sector reduces to a ferromagnetic Ising chain in a uniform 
transverse (pscudo)magnctic field: 

H T , cS = -J T £ t*t* +1 - h £ r* . (25) 

n n 

Notice that the restriction A <C J T , which was imposed 
in the derivation of the effective Hamiltonian in the spin 
sector, now can be released because the spin sector is 
assumed to be in the Neel phase. 

At h = J T , i.e. at A = A C 2 where A C 2 satisfies the 
equation 

Ac2 C(A c2 ) = J T , (26) 

the model (|25p undergoes a 2D Ising transitio n 27 i 31 to a 
massive disordered phase with (r^) = 0. This quantum 
critical point can be reached when A is further increased 
in the region A > A c x- It is clear from (|26p that X c2 is of 
the order of or greater than J T . It is reasonable to assume 
that for such values of A the Neel magnetization is close 
to its nominal value, £ ~ 1, implying that A C 2 ~ J T - We 
see that the two Ising transitions are well separated: 

A c2 /A cl ~ (J T /A S ) 1/2 » 1. (27) 

Thus, in the limit J T 3> A5, the ground-state phase 
diagram of the model (JTJ) consists of three gapped phases 



(a) (b) 



o 




FIG. 2: Schematic diagram of order parameters as functions 
of the SO coupling constant A. (a) Two Ising transitions in 
the J T 3> Ag limit, (b) A single Gaussian transition in the 
As S> Jt limit. These two scenarios correspond to path-1 
and path-2 in the phase diagram (Fig. 1), respectively. 



separated by two Ising criticalities, one in the spin sector 
(A = Ad) an d the other in the orbital sector (A = A C 2)- 
At < A < A c i the spin sector represents an anisotropic 
spin-liquid while in the orbital sector there is a Necl-likc 
ordering of the pseudospins: (— 1)™(t*) = 77^ (A) 7^ 0. 
At A c i < A < A C 2 the orbital degrees of freedom reveal 
their quantum nature: the onset of the spin Neel order 
(C 7^ 0) is accompanied by the emergence of the trans- 
verse component of the staggered pscudospin density: 
(— 1)"(t^) = r) x (X) 7^ 0. Upon increasing A, the stag- 
gered orbital order parameter r) undergoes a continuous 
rotation from the z-direction to ir-direction. At A = A c2 a 
quantum Ising transition takes place in the orbital sector 
where rf vanishes. At A > A C 2 both sectors are long- 
range ordered, with order parameters £, rf 7^ 0. The de- 
pendence of order parameters on A is schematically shown 
in Fig. 2(a); this picture is in full qualitative agreement 
with the results of the recent numerical studies^ 

The crossover between the small and large A limits 
studied in this section corresponds to path 1 on the phase 
diagram shown in Fig. [1] The path is located in the re- 
gion J T 3> A5. Starting from the massive phase I and 
moving along this path we first observe the spin-Ising 
transition (I — > II) to the Neel phase. Long-range or- 
dering of the spins induces quantum reconstruction of 
the initialy classical orbital sector (i.e. generation of a 
nonzero if). The orbital-Ising transition (II — > III) takes 
place inside the spin Neel phase. Of course, feedback ef- 
fects (that is, orbit affecting spin) become inreasingly im- 
portant upon deviating from the critical curve A5 J T ~ 1 
into phases II and III, especially in the vicinity of the 
orbital transition where the spin-orbit coupling is very 
strong, A ~ J T . In this region the behavior of the spin 
degrees of freedom is not expected to follow that of an 
isolated anisotropic spin-1 chain in the Neel phase since 
the effect of an "explicit" staggered magnetic field ~ Xif 
becomes important. We will see a pattern of such be- 
havior in the opposite limit of "heavy" spins, which is 
discussed in the next section. 
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IV. GAUSSIAN CRITIC ALITY AT J T <A S 

In this section we turn to the opposite limiting case: 
A 5 3> J T - Now the spin degrees of freedom constitute 
the "fast" subsystem and can be integrated out to gen- 
erate an effective action in the orbital sector. We will 
show that, in this regime, the intermediate massive phase 
where the orbital order parameter T] undergoes a continu- 
ous rotation from ry = (0, 0, rf) to rj = (rj x , 0, 0) no longer 
exists. Going along path 2, Fig.[T] which is located in the 
region As S> J T , we find that the two massive phases, I 
and III, are separated by a single Gaussian critical line 
characterized by central charge c = 1. On this line the 
vector rj vanishes, the orbital degrees of freedom become 
gapless and represent a spinless Tomonaga-Luttinger liq- 
uid characterized by power-law orbital correlations. 

At A = the spin-1 subsystem represents a disordered, 
isotropic spin liquid. Therefore the first nonzero correc- 
tion to the low-energy effective action in the orbital sector 
appears in the second order in A: 

ASW=-±(S St ) s (28) 
= 4 A "E / dT i / dr 2 (Sn(r 1 )S m (r 2 )) s TZ( n )T*(T 2 ), 

where (■ ■ -)s means averaging over the massive spin de- 
grees of freedom. According to the decomposition of the 
spin density, Eq. ((4]), the correlation function in ((29)) has 
the structure: 



(S,(r)S (0)) = (-l)7i(r/&) + / 3 (r/&) 



(29) 



Here £ s = v s /As is the spin correlation length and 
r = (v s t, x) is the Euclidian two-dimensional radius- 
vector, fi and f 2 are smooth functions with the following 
asymptotic behaviour— 

h(x) = dx-^e-*, f 2 (x) = C 2 x- l e- 2x (x » 1), (30) 

where C\ and C 2 are nonuniversal constants. DMRG 
calculations show^ that C 2 <§C C\\ for this reason the 
contribution of the smooth part of the spin correlation 
function can be neglected in (|29|) . 



Integrating over the relative time t_ = n — t 2 we find 
that the spin-orbit coupling generates a pseudospin xx- 
exchange with the following structure: 



h; = ^^(-i) ,+1 j;(!)tx +i 

n 1>1 



(31) 



Here the exchange couplings exponentially decay with 
the separation I, J' T {1) ~ (A 2 /Ag) exp(— Za /£s), so the 
summation in (|3ip actually extends up to I ~ £s/ao- In 
the Heisenberg model £s is of the order of a few lattice 
spacings, so for a qualitative understanding it would be 
sufficient to consider the I = 1 term as the leading one 
and treat the I = 2 term as a correction. Making a ir/2 



rotation in the pseudospin space, — »• t%, t% —> — t%, 
we pass to the conventional notations and write down the 
effective Hamiltonian for the orbital degrees of freedom 
as a perturbed XY spin-1/2 chain: 



H T — (JxT n T n+l + Jy T n T n+l) + ^t- 



where 



Hi. 



- V T X T X 

J x 'n 'n+2 



(32) 



(33) 



Here J y = J T , J x = J' T (l) > and J' x = J' r {2) > 0. By 
order of magnitude J' x < J x ~ A 2 / As- 

In the absence of the perturbation H' T , the model (|3"2")l 
represents a spin-1/2 XY chain which for any nonzero 
anisotropy in the basal plane (J x ^ J y ) has a Neel long- 
range order in the ground state and a massive excitation 
spectrum. This follows from the Jordan- Wigner trans- 
formation 

< = 24 a n - 1, r+ = r* +ir* - 2o+ e i7r ^ a >> (34) 

which maps the XY chain onto a model of complex spin- 
less fermions with a Cooper pairing^ 



H 



cir 



(Jx + Jy) ^2 ( a « an +! + h -°) 



+ (J x -J y )J2{ a i^+i+h.c). (35) 

n 

By increasing A (equivalently, decreasing J T ) the model 
(|35[) can be driven to a XX quantum critical point, J x = 



A 



A c ~ \/J T As, where the the system 



acquires a continuous U(l) symmetry. At this point the 
Jordan- Wigner fermions become massless and the system 
undergoes a continuous quantum transition. 

The transition is associated with reorientation of the 
pscudospins. Away from the Gaussian criticality the 
effective orbital Hamiltoian is invariant under Z 2 x Z 2 
transformations: t x — > —t x , — > — r*. In massive 
phases this symmetry is spontaneously broken. Making a 
back rotation from t v to t z we conclude that at J y > J x 
(A < A c ) ?f ^ 0, jf — 0, while at J y < J x (A > A c ) 
if = 0, rf 7^ 0. Both rf and rf vanish at the critical 
point, so contrary to the case J T 3> As, here there is no 
region of their coexistence. 



The passage to the continuum limit for the model ([32]) 
based on Abelian bosonization is discussed in Appendix 
[Bl There we show that the perturbation H' T adds a 
marginal four-fermion interaction g = J' x (2)/ttv <C 1 
to the free-fermion model (|B3[) . In the spin-chain lan- 
guage, this is equivalent to adding a weak ferromag- 
netic zz-coupling. In the limit of weak XY anisotropy, 
| A — A c |/A c <C 1, the low-energy properties of the orbital 
sector are described by a quantum sine-Gordon model 
(all notations are explained in Appendix [Bjl 



*-5 



atj 2 + 1 (a**) 



— cos v^e, (36) 

ITCX. 
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where 



7 ~ Jt 



A - Ac 

A, 



K =l + 2.g + 0(<7 2 )- 



(37) 



The U(l) criticality is reached at A = A c where, due to a 
finite value of g, the orbital degrees of freedom represent 
a Tomonaga-Luttinger liquid. Close to the criticality, the 
spectral gap in the orbital sector scales as the renormal- 
ized mass of the sine-Gordon model (I3"6l) : 



A - A, 



Ac 



(38) 



Strongly fluctuating physical fields acquire coupling 
dependent scaling dimensions. In particular, according to 
the bosonization rules^ the staggered pseudospin den- 
sities are expressed in terms of the vertex operators, 



(-!)"< 



n x {x) 
n z (x) 



sin v / 7rO(x) 
cos ^/tt@(x) 



(39) 



both with scaling dimension d = l/AK. This anomalous 
dimension determines the power-law behaviour of the av- 
erage staggered densities close to the criticality: 



rf(A) 
rf(A) 



A < A c 
A > A, 



(40) 



A finite staggered pseudospin magnetization rf at 
A > A c generates an effective external staggered mag- 
netic field in the spin sector: 



H = Hq 



H' s , 



H's 



-^£(-1)"^, ( 41 ) 



where hs = —\ri x . The spectrum of the Hamiltonian H 
is always massive. This can be easily understood within 
the Majorana model ((2|). According to ([6]), in the con- 
tinuum limit, the sign-alternating component of the spin 
magnetization, N 3 ~ (—1)™^, can be expressed in terms 
of the order and disorder fields of the degenerate triplet 
of 2D disordered Ising models: N 3 ~ fii^os. In the 
leading order, the magnetic interaction H' s gives rise to 
an effective magnetic field h% = hs(/J>xfJ>2) applied to the 
third Ising system: /13O3. The latter always stays off- 
critical. 

Since in the Haldane phase the spin correlations are 
short-ranged, close to the transition point the induced 
staggered magnetization £ can be estimated using linear 
response theory. Therefore, at < A — A c A c , ( follows 
the same power-law increase as that of rf but with a 
smaller amplitude: 



C hi. ( Jt 

4 ~ A s ~ \A S 



1/2 



A — A c 
A c 



1/4K 



(42) 



So, in the part of the phase C, Fig. [TJ where A 5 ^> J T , 
the ^-orbital order, being the result of a spontaneous 



breakdown of a Z2 symmetry t x — > — r%, acts as an ef- 
fective staggered magnetic field applied to the spins and 
induces their Neel alignment. This fact is reflected in a 
coupling dependent, nonuniversal exponent 1/AK char- 
acterizing the increase of the staggered magnetization at 
A > A c . The order parameters as functions of A in the 
As ^> J T limit is schematically shown in Fig. 2(b). 

As already mentioned, the absence of a small parame- 
ter in the regime of strong hybridization, J T ~ J$ ~ A, 
makes the analysis of the phase diagram in this region not 
easily accessible by analytical tools. Nevertheless some 
plausible arguments can be put forward to comment on 
the topology of the phase diagram. It is tempting to 
treat the curve J T Ag/A 2 ~ 1 as a single critical line go- 
ing throughout the whole phase plane (J T /X, Ag/A). If 
so, we then can expect that there exists a special sin- 
gular point located in the region J T Ag/A 2 ~ 1. This 
expectation is based on the fact that at J T S> A s limit 
the transition is of the Ising type and the spontaneous 
spin magnetization below the critical curve follows the 
law £ ~ (A — Ad) 1 / 8 with a universal critical exponent, 
whereas at J T -C A s the spin magnetization has a differ- 
ent, nonuniversal exponent, £ ~ (A — Ac) 1 ^ 4 ^. Continu- 
ity considerations make it very appealing to suggest that 
at the special point the Tomonaga-Luttinger liquid pa- 
rameter takes the value K — 2, and the two power laws 
match. Since the central charges of two Ising and one 
Gaussian criticalitics satisfy the relation 1/2 + 1/2=1, 
the singular point must be a point where the two Ising 
critical curves merge into a single Gaussian one. 



V. DYNAMICAL SIN SUSCEPTIBILITY AND 
NMR RELAXATION RATE IN THE VICINITY 
OF GAUSSIAN CRITICALITY 



It may seem at the first sight that, in the regime 
As 3> J, the spin degrees of freedom which have been 
integrated out remain massive across the orbital Gaus- 
sian transition, and the spectral weight of the staggered 
spin fluctuations is only nonzero in the high-energy re- 
gion ui ~ As- However, this conclusion is only correct 
for the zeroth-order definition of the spin field No (a;), 
given by Eq. with respect to the spin-orbit inter- 

action. In fact, the staggered magnetization hybridizes 
with low-energy orbital modes via SO coupling already in 
the first order in A and thus acquires a low-energy pro- 
jection which contributes to a nonzero spectral weight 
displayed by the dynamical spin susceptibility at ener- 
gies well below the Haldane gap. 

To find the low-energy projection of the field N z (r), we 
must fuse the local operator N§ (r) with the perturbative 
part of the total action. Keeping in mind that close to 
and at the Gaussian criticality most strongly fluctuating 
fields are the staggered components of the orbital po- 
larization, we approximate the SO part of the Euclidian 
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action by the expression 

S St * — / d 2 r N z (r)n*(v), (43) 

where r = (vst,x) is the two-dimensional radius vector 
(here r is the imaginary time). We thus construct 

N 7 p{r) = (e- Ss ^N z (r)) 

= N§(v) - ^° [ d 2 n(^ WS(*i))s » x (ri) 

+ 0(A 2 ), (44) 

where averaging is done over the unperturbed, high- 
energy spin modes. For simplicity, here we neglect the 
anisotropy of the spin-liquid phase of the S=l chain and 
use formula ([8]). The spin correlation function is short- 
ranged. Treating the spin correlation length £s ~ ug/Ag 
as a new lattice constant (new ultraviolet cutoff) and be- 
ing interested in the infrared asymptotics |r| » £g, we 
can replace in (|4"4")l n x (rx) by n x (r). The integral 

J d 2 p (NZ(p)N%(0)) s (45) 

~ -Nts) m / dp p^fp e-oK* ~ (Cs/ao) 5/4 . 

So the first-order low-energy projection of the staggered 
magnetization is proportional to 

N z P (r) ~ (-) 1 n*(r). (46) 
A s V a o/ 

This result clarifies the essence of the hybridization effect: 
close to the Gaussian criticality the spin fluctuations ac- 
quire a finite spectral weight in the low-energy region, 
lj <C Ag, q ~ 7r, which is contributed by orbital fluctu- 
ations and can be probed in magnetic inelastic neutron 
scattering experiments and NMR measurements. 

Away from but close to the Gaussian criticality the 
behavior of the dynamical spin susceptibility Qmx(q, uS) 
is determined by the excitation spectrum of the sine- 
Gordon model for the dual field, Eq. ([3l)|) . Since K > 1, it 
consists of kinks, antikinks carrying the mass M or b, and 
their bound states (breathers) with masses (see e.g. Ref. 

S3) 

Mj = 2M orb sin(7rj/2^), 

j = 1,2,... v - 1, v = 2K-l (47) 

Since K — 1 = 2g is small, there will be only the first 
breather in the spectrum, with mass Mi = 2A/ orD (l — 
2ir 2 g 2 ). The sine-Gordon model is intcgrable, and the 
asymptotics of its correlation functions in the massive 
regime have been calculated using the form-factor ap- 
proach (sec for a recent review |35|). Here we utilize 
some of the known results. At A < A c the operator 
n x ~ sin y^O has a nonzero matrix element between 



the vacuum and the first breather state. This form- 
factor contributes to a coherent peak in the dynamical 
spin susceptibility at frequencies much smaller than than 
the Haldane gap: 

3m X (g, W ,T = 0) = A(X/ A s ) 2 S[uj 2 — (q — ir) 2 v 2 — Alf] 
+ 3TOXco„t(<7,w,T = 0). (48) 

Here A is a constant and the second term is the contribu- 
tion of a multi-kink continuum of states with a threshold 
at u> = 2M or b- At A > A c the spectral properties of 
the operator cos y/irO coincide with those of the opera- 
tor sm^/nQ at A < A c . For symmetry reasons^, this 
operator does not couple to the first breather, so that 
at A > A c 3mx(?, w) will only display the kink-antikink 
scattering continuum. 

We see that, due to spin-orbit hybridization effects, the 
spin sector of our model loses the properties of a spin liq- 
uid already in a noncritical orbital regime. This tendency 
gets strongly enhanced at the orbital Gaussian criticality 
(Morb — > 0) where all multi-particle processes merge, and 
the spin correlation function exhibits an algebraically de- 
caying asymptotics 

(N*(r)N*(0)) ~ (iV* (r)JVJ,(0)) ~ (£) 

(49) 

implying that the spin sector of the model becomes remi- 
niscent of Tomonaga-Luttinger liquid. In this limit (here 
for simplicity we consider the T = case) the dynamical 
spin susceptibility is given by^l 

Zm X (q, u, T = 0) ~ (A/A s ) 2 [u 2 - v 2 (q - tt) 2 ] . 

(50) 

The NMR relaxation rate probes the spectrum of local 
spin fluctuations 

±- = A 2 T lim - V V z (?, w, T) 

1 1 oj->0 Ul z — ' 

1 

where A is an effective hypcrfine constant. In spin- 
liquid regime of an isolated spin-1 chain, the existence of 
a Haldane gap makes 1/T\ exponentially suppressed^: 
1/T\ ~ exp(— 2As/T). The admixture of low-energy or- 
bital states in the spin-fluctuation spectrum drastically 
changed this result. A simple power counting argument*^ 
leads to a power-law temperature dependence of the 
NMR relaxation rate: 

k~ A2 (€) T ^ (5i) 

This result is valid not only exactly at the Gaussian crit- 
icality but also in its vicinity provided that the tempera- 
ture is larger than the orbital mass gap. By construction 
(see the preceding section) K > 1. This means that the 
exponent 1/2 A" — 1 is negative and the NMR relaxation 
rate increases on lowering the temperature. It is worth 
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noticing that such regimes are not unusual for Tomonaga- 
Luttinger phases of frustrated spin-1/2 ladders^ For our 
model, such behavior of 1/Ti would be a strong indica- 
tion of an extremely quantum nature of the collective 
orbital excitations^ 



VI. BEHAVIOR IN A MAGNETIC FIELD: 
QUANTUM ISING TRANSITION IN ORBITAL 
SECTOR 

We have seen in Sec lIIII that, due to spin-orbit cou- 
pling, the Neel ordering of the spins is accompanied by 
the emergence of quantum effects in the orbital sector: 
the classical orbital Ising chain transforms to a quan- 
tum one. In this section we briefly comment on a similar 
situation that can arise upon application of a uniform 
external magnetic held h. 

Since the spin-1 chain is massive, it will acquire a finite 
ground-state magnetization (S z ) only when the magnetic 
held, h, is higher than the critical value h c \ ~ Ag, corre- 
sponding to the commensurate-incommensurate (C-IC) 
transition. According to the definition (|5|), a uniform 
magnetic field along the z-axis, H ma ,g = —hl z , mixes 
up a pair of Majorana fields, £ x and £ 2 , and splits the 
spectrum of S z = ±1 excitations (the S z = modes are 
unaffected by the field). At h = h c i the gap in the spec- 
trum of the S z = 1 excitations closes, and at ft, > h c i 
these modes condense giving rise to a finite magnetiza- 
tion. Once (S z ) 7^ 0, the effective Hamiltonian of the 
T-chain becomes 

H T = J T ^rX z +1 -A T ^<, A T = X(S Z ). (52) 

n n 

Here we ignore the fluctuation term that couples t x to 
AS Z — S z — \S%)- 

One should keep in mind that there exists the sec- 
ond C-IC transition at a higher field h C 2 associated with 
full polarization of the spin-1 chain. To simplify further 
analysis, let us assume that the range of magnetic fields 
h c i < h < h C 2, where an isolated spin-1 chain has an in- 
commensurate, gapless ground state, is sufficiently broad. 
This can be easily achieved in the biquadratic model ^ 
with /3 1, in which case the Haldane gap - and hence 
h c i ~ is small, and the effects associated with the second 
C-IC transition can be neglected. 

Now, by increasing the magnetic field h in the region 
h > h c i, the effective orbital chain (|52p can be driven 
to an Ising criticality. The induced transverse "magnetic 
field" A T is proportional to a nonzero magnetization of 
the spin-1 chain. If A/ J T is large enough, then upon in- 
creasing the field the effective quantum Ising chain ([5^)1 
can reach the point A T (h*) = J T where the Ising transi- 
tion occurs. This will happen at some field h = h* > h c i- 
In the region \h — h*\/h* <C 1 the quantum Ising T-chain 
will be slightly off-critical. Due to the SO coupling, these 



massive orbital excitations will interact with the gap- 
less S z = ±1 spin modes. However, this interaction can 
only give rise to the orbital mass renormalization (i.e. a 
small shift of the Ising critical point) and a group veloc- 
ity renormalization of the spin-doublet modes. For this 
reason we do not expect the aforementioned spin-orbital 
fluctuation term to cause any qualitative changes. 

The above discussion reveals an interesting fact: a suf- 
ficiently strong magnetic field acting on the spin degrees 
of freedom can affect the orbital structure of the chain 
and drive it to a quantum Ising transition. The difference 
with the situation discussed in Sec lIIIl is that the external 
magnetic field induces a uniform spin polarization which, 
in turn, gives rise to a uniform transverse orbital order- 
ing (t%) 7^ 0. Thus, the classical long-range orbital order 
( r n) = (~l) n? 7 z 5 present at h < h* , disappears in the 
region h > h* , where the orbital degrees of freedom are 
characterized by a transverse ferromagnetic polarization, 
(r x ) ± 0. 

VII. CONCLUSION AND DISCUSSION 

In this paper, we have proposed and analyzed a ID 
spin-orbital model in which a spin-1 Haldane chain is lo- 
cally coupled to an orbital Ising chain by an on-site term 
Xt x S z originating from rclativistic spin-orbit (SO) in- 
teraction. The SO term not only introduces anisotropy 
to the spin sector, but also gives quantum dynamics to 
the orbital degrees of freedom. We approach this prob- 
lem from well defined limits where either the spin or the 
orbital sector is strongly gapped and becomes a 'fast' 
subsystem which can be integrated out. By analyzing 
the resultant effective action of the remaining 'slow' de- 
grees of freedom, we have identified the stable massive 
and critical phases of the model which are summarized 
in a schematic phase diagram shown in Fig. 1. 

In the limit dominated by a large orbital gap, i.e. J T 3> 
As, integrating out the orbital variables gives rise to an 
easy-axis spin anisotropy D(S Z ) 2 where D ~ — A 2 /J T . 
As A increases, the disordered Haldane spin liquid un- 
dergoes an Ising transition into a magnetically ordered 
Neel state. The presence of antifcrromagnctic spin order 
( in the Neel phase in turn generates an effective trans- 
verse field h ~ X( acting on the orbital Ising variables. 
The orbital sector which is described by the Hamiltonian 
of a quantum Ising chain reaches criticality when h = J T . 
In between the two Ising critical points lies an interme- 
diate phase (phase II in Fig. 1) where both Ising order 
parameters i] x and rf are nonzero. Such a two-stage or- 
dering scenario illustrated by path 1 in the phase diagram 
(Fig. 1) has been confirmed numerically by recent DMRG 
calculations^ Interestingly, the orbital Ising transition 
can also be induced by applying a magnetic field to the 
spin sector. As the field strength is greater than the Hal- 
dane gap, a field-induced magnon condensation results 
in a finite magnetization density (S z ) in the linear chain. 
Thanks to the SO coupling, the orbital sector again ac- 
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quires a transverse field h ~ A(S' Z ) and becomes critical 
when h = J T . 

A distinct scenario of the orbital reorientation tran- 
sition rf — > rf occurs in the opposite limit A5 3> J T . 
This time we integrate out the fast spin subsystem and 
obtain a perturbed spin- 1/2 XY Hamiltonian for the or- 
bital sector. The effective exchange constants are given 
by J x ~ A 2 /As and J y — J T . As A is varied, the orbital 
sector reaches a Gaussian critical point when J x = J y , 
at which the system acquires an emergent U(l) sym- 
metry. The orbital order parameter goes directly from 
r\ = (0, 0, ij z ) to {rf , 0, 0) in this single-transition scenario 
(illustrated by path 2 in Fig. 1). Both order parameters 
rf and rf vanish at the critical point. We have shown 
that spin-orbital hybridization effects near the Gaussian 
transition lead to the appearance of a non-zero spectral 
weight of the staggered spin density well below the Hal- 
dane gap - the effect which can be detected by inelastic 
neutron scattering experiments and NMR measurements. 

The stability analysis of the orbital Gaussian criti- 
cality in the original lattice model ([T]), done in Ap- 
pendix |Bl has shown that this critical regime is pro- 
tected by the t z — ¥ —t z symmetry of the underlying 
microscopic model. This symmetry will be broken in 
the presence of an orbital field 5 T n which removes 
degeneracy between the local orbitals d zx and d yz and 
adds a "magnetic" field along the y-axis in the effective 
XY model ([52")) . Such perturbation will drive the orbital 
sector away from the Gaussian criticality. The same 
argument applies to a perturbation with the structure 
P T n ^n T n which a i so breaks the aforementioned sym- 
metry. Integrating over the spins will generate an extra 
term ~ A/3 J2n ( T n T n+i + T n T n+i) which, in the contin- 
uum limit, translates to A/3 sin a/47t6. As explained in 
Appendix [B] such perturbation will keep the orbital sec- 
tor gapped with coexisting rj x and rf orderings. 

Since the analysis presented in this paper is done in the 
limiting cases, precise predictions on the detailed shape 
of the phase diagram or on the behavior of correlation 
functions in the regime of strong hybridization of spin 
and orbital degrees of freedom, where all interactions in- 
cluded in the model are of the same order, are beyond our 
reach and require further numerical calculations. On the 
other hand, the continuity and scaling analysis allow us 
to believe that the global topology of the phase diagram 
and character of critical lines are given correctly. Finally 
the spin-orbital model Eq. ([1]) can be generalized to the 
zigzag geometrical where two parallel spin-1 chains are 
coupled to a zigzag Ising orbital chain via on-site SO in- 
teraction. The zigzag case is closely related to the quasi- 
1D compound CaV204. While the two-Ising-transitions 
scenario is expected to hold in the J T ^> A5 regime, the 
counterpart of Gaussian criticality in the zigzag chain 
remains to be explored and will be left for future study. 
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Appendix A: Ising correlation function 

In this Appendix we estimate the correlation func- 
tion r^(r) = (t^(t)t^(O)), where the averaging is per- 
formed over the ground state of the Ising Hamiltonian 

It proves useful to make a duality transformation: 

T Z T Z — ll X T X — ll Z ll Z 

The new set of Pauli matrices fx^ represents disorder op- 
erators. The Hamiltonian and correlation function be- 
come: 

H^J T J2Vn, (Al) 

n 

K Z m (r) -> (^(tK+i(t)/4(0)^ +1 (0)>. (A2) 

The most important fact about the dual representation 
is the additive, single-spin structure of the Hamiltonian: 
the latter describes noninteracting spins in an external 
"magnetic field" J T . Notice that by symmetry (/x z ) = 0. 
Therefore the correlation function in (|A2[) has an ultralo- 
cal structure: 

KM = S nm Y 2 (r), Y(t) = «(r)<(0)). (A3) 

The time-dependence of the disorder operator can be ex- 
plicitly computed, 

f£( r ) = e Tj T^n z n e- Tj -r^ = /4cosh(2J T T)-i/^ sinh(2 J t t) 

Therefore (below we assume that r > r') 

F(r-r') = cos1i2J t (t-t') + (/z x )sinh 2 J t (t-t') 
= exp[-2|J T |(r-r')]. (A4) 

Here we used the fact that, in the ground state the Hamil- 
tonian H T , (fi x ) = — sgn J T . Thus, as expected for the 
ID Ising model, the correlation function r^(r) is local 
in real space and decays exponentially with r: 

r^(r) = <5 ram exp(-4J ± |r|). (A5) 
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Appendix B: Perturbed XY chain, Eq. (f32|) 



In this Appendix we analyze the perturbation (|33[) to 
the XY spin chain (|3"2"]l and show that at the XX point 
it represents a marginal perturbation which transforms 
the free-fermion regime to a Gaussian criticality describ- 
ing a Luttinger-liquid behavior of the orbital degrees of 
freedom. 



Using the Jordan- Wigner transformation (|34[) we 
rewrite (J33J) as H' T = H[+H' 2 , where 

H , = JL^j2(aia n+2 + h.c.)(ai + ian + i-h(Bl) 



H' = 



4(2) 



(44+2 + h.c.)(a f n+1 a n+1 - i). (B2) 



Assuming that | J x — J y | , J' x <C J x + J y , we pass to a con- 
tinuum description of the XY chain in terms of chiral, 
right (R) and left (L), fermionic fields based on the de- 
composition (to simplify notations we set here ag = 1): 
a n —> (— i) n R(x) + i n L(x). Then the Hamiltonian density 
of the XY model takes the form: 



-Hxy(x) = -iv (R)d x R - tfd x L) 



2i 7 (R)]J - h.c.) , 

(B3) 



where 



Standard rules of Abelian 



1 = Jx »y 

bosonization^i transform (|B3|) to a quantum sine-Gordon 
model: 



T-Lxy{x) 



n- 



20 



cosV47r8, (B4) 



where v = 2(J X + J y )ao is the Fermi velocity, Tl(x) — 
d x Q(x) is the momentum conjugate to the scalar field 
= $r(x) + <$> L (x), and &(x) = -$r(x) + $ L (x) 
is the field dual to Here $>r.l(x) are chiral 

components of the scalar field. Using the fact that 
the fermions are spinless, one can impose the condi- 
tion $l(cc')] = i/4 and thus make sure that the 
bosonization rules correctly reproduce the anticommuta- 
tion relations {R(x),L(x')} = {R(x), L^{x')} = 0. An 
explicit introduction of the so-called Klein factors be- 
comes necessary when bosonizing fermions with an inter- 
nal degree of freedom, such as spin 1/2, chain index etc, 
which is not the case here. 

Let is find the structure of the perturbation (|33[) in the 
continuum limit. First of all we notice that 



1 n+l a n+l 



1/2 



(: R^R: + : L ] L :) + (-l) n+1 (R^L + L ] R) 



(-1)" 



sin v 47r$ 



(B5) 



Similarly 

4 a «+2 + h.c. 

->■ -2 [(: R^R: + : L^L :) + (-l) n (R?L + L^R)] 



= -2 



1 f— 1)™ 

—=d x <S> - — sin V47r$ 

■J 7r 7ra 



(B6) 



Dropping Umklapp processes R^{x)R^{x + a)L(x + 
a)L(x) + h.c. ~ cos vT67r<f> as strongly irrelevant (with 
scaling dimension 4) at the XX criticality and ignoring 
interaction of the fermions in the vicinity of the same 
Fermi point, we find that 

(4«n+2 + h.c.)(al +1 a n+1 - 1/2) 

smooth 

-S- -8 : R^R :: L^L := 2 [IT 2 - (d x $) 2 ] . (B7) 

We see that the perturbation H[ generates a marginal 
four-fermion interaction to the free-fermion model (|B3[) . 
thus transforming the model (|32|) to an XYZ model with 
a weak ferromagnetic (zz)-coupling. This interaction can 
be incorporated into the Gaussian part of the bosonic 
theory (|B4I) by changing the compactification radius of 
the field $: 



H 



Hxy + n'i 

KU 2 + 



— cosVA^Q. (B8) 



Here u is the renormalizcd velocity and K is the inter- 
action constant which at J' x <C (J x + J y ) is given by 
K = 1 + 2g + 0(g 2 ), where g = J x (2)a /nv < 1. 

Now we turn to H' . We have: 



at a; 



h. 



-T ( B9 ) 

a; + a ) + L\x)R' i {x + a) + h.c] 
+(-1)" [-Rt^-Rtfa; + a) + L^L^x + a) + /i.e.] . 



Bosonizing the smooth term in the r.h.s. of (|B10[) one 
obtains d x <& cos \/47r0. Bosonizing the staggered term 
yields sin \/Att<& cos y/AirQ. Using the OPE 



sin VAtt^(x) sin V / 47r$(x + a) 



= const - na 2 (d x §f 



- cos Vl67r(f>, 
2 



we find that, in the continuum limit, the Hamiltonian 
density R! 2 is contributed by the operators cos a/Itt© and 
(d x <&) 2 cos %/4~7r0 (as before, we drop corrections related 
to Umklapp processes). The former leads to a small ad- 
ditive renormalization of the fermionic mass 7 and thus 
produces a shift of the critical point. The latter repre- 
sents an irrelevant perturbation (with scaling dimension 
3) at the XX criticality. In a noncritical regime it renor- 
malizes the mass and four-fermion coupling constant g. 

Considering the structure of the remaining terms in 
the expansion pip one arrives at similar conclusions. 
Here a remark is in order. The only dangerous perturba- 
tion which would dramatically affect the above picture 
is sin-\/47r0. The presence of two nonlinear terms in 
the Hamiltonian, 7 cos VAttQ + S sin y/AnQ, would make 
the fermionic mass equal to \/(A — A c ) 2 + S 2 . The Gaus- 
sian criticality in this case would never be reached, the 
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model would always remain massive, and nonzero stag- 
gered pseudospin densities, rf and rf , would coexist in 
the whole parameter range of the model. 

Fortunately, the appearance of the operator sin V^r© 
is forbidden by symmetry. The initial Hamiltonian ([T]) 
is invariant under global pseudospin inversion In the z- 
component only: — > —t*. After rotation t z — > t v 



this translates to — > — t*. Using the bosonized ex- 
pressions ([3"9")) for the staggered pseudospin densities we 
find that the corresponding transformation of the dual 
field is 6 — > ^/tt — and so the bosonized Hamiltonian 
density must be invariant under this transfomation. This 
explains why the operator sin y/4nQ cannot appear in the 
effective continuum theory. 
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